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Pulse-splitting in light propagation through A^-type atomic media 
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We investigate the spatio-temporal evolution of a Gaussian probe pulse propagating through a 
four-level Wtype atomic medium. At two-photon resonance of probe-and control fields, weaker 
probe pulses may propagate through the medium with low absorption and pulse shape distortion. 

In contrast, we find that increasing the probe pulse intensity leads to a splitting of the initially 
Gaussian pulse into a sequence of subpulses in the time domain. The number of subpulses arising 
throughout the propagation can be controlled via a suitable choice of the probe and control field 
parameters. Employing a simple theoretical model for the nonlinear pulse propagation, we conclude 
that the splitting occurs due to an interplay of Kerr nonlinearity and group velocity dispersion. 

PACS numbers: 42.50.Gy, 32.80.Qk, 42.65.-k 


I. INTRODUCTION 

Atomic coherence and quantum interference effects are 
key drivers for the control of light-matter interactions. 
Among the most prominent effects are electromagneti- 
cally induced transparency (EIT) [BIS] and related phe¬ 
nomena like the slowing and stopping of light [aHl!, the 
coherent storage and retrieval of light HHn], and sta¬ 
tionary light [Mil]- EIT can also be used to enhance 
the nonlinear susceptibility [B EMSS]. However, con¬ 
ventional EIT schemes are limited in their enhancement 
of nonlinear effects due to the presence of considerable 
absorption [BUI]- Multilevel EIT-related atomic sys¬ 
tems are a possible route to circumvent these limita¬ 
tions. In this spirit, Schmidt and Imamoglu proposed 
a four-level TV-configuration system in which the stan¬ 
dard EIT level scheme is augmented by an additional 
off-resonantly-driven transition for the enhancement of 
the Kerr-nonlinearity |26j . This Kerr-nonlinearity was 
experimentally observed in m and has many interesting 
applications. For example, the A^-type level system has 
been shown to exhibit cross-Kerr interactions [28] as well 
as the phenomena of optical bistability |53| and optical 
switching [30] in cavity systems and hollow-core optical 
fibres m- On the quantum level, the iV-level system 
gives rise to nonlinear optical interactions that can be 
conservative |3S] or dissipative [33] • These two-particle 
interactions were employed for the realization of strongly 
correlated polariton systems [32H34] . 

From a practical point of view, it is important to mea¬ 
sure the effect of the Kerr-nonlinearity on the probe held. 
An obvious and common way of probing an optical sys- 
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tern is by studying its response to an optical pulse. In 
particular, the polariton systems in [331 - 133] involve a se¬ 
quence where an input pulse is stored, processed and then 
released for readout. An evident question is thus how 
the Kerr-nonlinearity in the A^-level system modihes the 
propagation of a weak, classical probe pulse. However, 
up to now this problem has not received much atten¬ 
tion [331I37| - 

Motivated by this, here, we investigate the pulse prop¬ 
agation through a coherent medium in A^-type conhgu- 
ration, see Fig. [^ We solve the Maxwell-Bloch equa¬ 
tions governing the propagation dynamics numerically 
and consider resonant driving of the Raman transition 
within the A sub-system. As a hrst result, we hnd 
that a weak probe may propagate through the medium 
with low absorption and pulse shape distortion, provided 
that the Kerr detuning of the probe field is sufficiently 
large. However, at stronger probe pulse intensities, the 
numerical simulations predict a splitting of the initially 
Gaussian pulse into sub-pulses, which subsequently prop¬ 
agate through the medium without further splittings. 
To interpret these results, we derive a simple analyti¬ 
cal model for nonlinear pulse propagation through the 
A^-type medium. From this model, we conclude that the 
pulse splitting is a direct manifestation of the self-phase 
modulation associated with the Kerr nonlinearity in the 
presence of group velocity dispersion. The group velocity 
dispersion can be adjusted via the detuning of the probe 
and control helds with respect to the excited state of the 
A sub-system. Measuring the temporal pulse shape of 
the probe pulse at the output thus indicates the strength 
of the Kerr nonlinearity experienced inside the medium. 
A similar pulse splitting phenomenon should be observ¬ 
able in other systems with EIT-enhanced nonlinearities 
giving rise to self-phase modulation effects. For example, 
EIT in Rydberg systems [33] can give rise to giant pho¬ 
ton non-linearities mediated by the interaction between 
the Rydberg atoms [dOHiSj . 

This paper is organized as follows. In Sec. [H] we intro- 
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FIG. 1: (Color online) (a) A block diagram of the system of 
interest. We assume that the control field intensity is constant 
along the propagation direction, (b) Schematic representation 
of four-level atoms in A^-type configuration. We consider two 
J = 1/2 Zeeman manifolds with excited states |1) = |me = 

— 1/2) and |2) = |me = 1/2) and ground states |3) = \mg = 

— 1/2) and |4) = \mg = 1/2). The CT“-polarized control field 
Ec couples to the transition |1) -H- |4), and the rr-polarized 
probe field Ep interacts with the |1) —>■ |3) and |2) —>■ |4) 
transitions. 

duce the Maxwell-Bloch equations governing the pulse 
propagation in the considered atomic four-level medium. 
In Sec. m the propagation dynamics of the pulse is stud¬ 
ied numerically at different parameter conditions. In or¬ 
der to find a qualitative interpretation of our numerical 
findings, we introduce a model for the description of pulse 
propagation in nonlinear media in Sec. |IV| Section |V] pro¬ 
vides a summary and discussion of our results. 

II. THEORETICAL MODEL 

We consider a homogeneous atomic medium interact¬ 
ing with two laser fields as shown in Fig. [^a). The weak 
probe and the strong control fields propagate in perpen¬ 
dicular directions. The atomic four-level system in N- 
type configuration is depicted in Fig. [^b). In the fol¬ 
lowing, we assume that this level scheme is realized with 
two J = 1/2 Zeeman manifolds via polarization selection. 
The weak probe field with central frequency ujp couples to 
the TT-transitions [|1) O |3) and |2) o |4)] and is defined 
as 

Ep{t) = + c.c., (1) 

where e,^. is the polarization unit vector, kp = uip/c and 
c is the speed of light. The strong CT“-polarized control 
field with frequency oJc couples to the |1) O |4) transition 
and can be written as 

£;c(t)=e,-4e*('==^-'"=*)+c.c., (2) 

where = uJc/c. In electric-dipole and rotating-wave 
approximation, the semiclassical Hamiltonian of the sys¬ 
tem is 

= -;i(Ap-Ac)|4)(4|-;iAp|l)(l| 

— li(Ap + Ap/— Ac)|2)(2| 

- [(|2)(4|-|l)(3|)Hp + |l)(4|H,+ H.C.] . (3) 


The Rabi frequencies of the probe and control fields in 
Eq. (§ are defined as 

np = ^^^£p , (4a) 

O, = , (4b) 

where d.y is the electric dipole moment for the transition 
between levels |i) and \j). Note that ^24 = —dis accord¬ 
ing to the Clebsch-Gordan coefficients of the considered 
J=l/2e-)'J = l/2 level scheme. The detunings of 
probe and control fields are given by 

Ap = ujp- wia , (5a) 

Ap/ =ujp- UJ24 , ( 5 b) 

Ac = <jJc- Wi4 ■ (5c) 


We model the quantum dynamics of the atoms using the 
master equation approach. The matrix elements of the 
atomic density operator p obey the following equations. 


Pll 

= 

-(7p + 7s)pii - 'i^pP3i 

+ i^*pPl3 




+ iHcP4i — i£l*pi4 , 


(6a) 

P 12 

= 

-[7p + 7s - *(Ac - Ap- 

-)]/0l2 + 

P42 



— iflpP32 — i£lpPl4, 


(6b) 

Pl3 

= 

-[(7p + 7s)/2 - iAp]pii 

j + i£lcP43 




— iflpipss — Pii), 


(6c) 

Pl4 

= 

-[(7p + 7s)/2 - iAc]pu 

t — i£lp{p34 

+ P 12 ) 



+ 1Hc(p44 — Pll), 


(6d) 

P 22 

= 

-(7p + 7s)P22 + i£ipP42 

- i%P24 , 

(6e) 

P 23 

= 

-[(7p + 7s)/2 + *(Ac- 

■ (Ap -b Ap 

’'))] P 23 



-|- iflpP43 + iflpP2i , 


(6f) 

P 24 

= 

- [(7p + 7s)/2 - *Ap/] P24 + i^pip44 - P 22 ) 



— i£lcP21 , 


(6g) 

P33 

= 

7pPll + 7sP22 - i^pPl3 - 

f i£tpP3i , 

(6h) 

P34: 

= 

[F -1- i{Ap — Ac)] P34 — 

i£l*pPi4 




i£tpP32 


(6i) 

where 

7p 

[7s] is the decay rate 

of the TT 

[cr] transi- 


tion. These decay constants are related through the 
Clebsch-Gordan coefficients of the level scheme and obey 
7s = 27p = 27/3, where 7 is the total decay rate of each 
excited state. The decay rate of the ground state coher¬ 
ence is denoted by F. The medium polarization induced 
by the probe field can be expressed in terms of the off- 
diagonal density matrix elements pis and p24 as 

P = Af + d42P2A) 1 (7) 

where Af is the number density of the medium. The 
spatio-temporal evolution of the probe pulse through the 
medium is governed by Maxwell’s equations. In the 
slowly varying envelope approximation, we find 

/ d 1 d\ 
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FIG. 2: (Color online) Intensity of a weak probe 

pulse as a function of position and time for different detun¬ 
ing conditions. The parameters in (a) and (b) are Q.c = 7 , 
fip = Vo.0157, CTp = 30/7, and r = 0. 



where 




ry = 7 


Stt 


(9) 


is the coupling constant and A is the wavelength of the 
transition |1) o |3). 

In order to facilitate the numerical integration of the 
coupled Maxwell-Bloch equations in Eqs. § and (|^, we 
introduce a co-moving coordinate system 


T = t-, ( = X. (10) 

c 


This change of coordinates reduces the partial differen¬ 
tial equation in Eq. (|^ to an ordinary differential equa¬ 
tion with respect to the independent variable V 


df2p 

— = IT ] (p 24 - P13) 


( 11 ) 


Typical parameter values are Af = 5x lO^^atoms/cm^, 
A = 780 nm, and 7 = 27r x 10® Hz. Throughout this 
paper we assume that the depletion of the strong control 
field throughout its propagation can be neglected such 



FIG. 3: (Color online) Probe pulse intensities at the medium 
exit at 7C/7 = 200. Parameters are as in Fig. To facil¬ 
itate the analysis of the pulse widths, all pulses are shown 
normalized to their respective peak values. 


that its intensity is independent of position. We further 
choose values of Ap and Ac compatible with the reso¬ 
nant Raman condition Ap = Ac, and hence we define 

A = Ap = Ac ■ (12) 


III. NUMERICAL ANALYSIS 


We start our discussion with a numerical analysis of 
the probe field propagation through the atomic four-level 
medium. To this end, we assume that the shape of the 
input pulse at the medium boundary at C = 0 is given by 

Hp(C = 0,r) = HOe , (13) 

where and Up are the amplitude and temporal width of 
the Gaussian pulse, respectively. Furthermore, all atoms 
are assumed to be in the ground state |3) initially. 

Our choice of parameters is guided by the following 
considerations. For suitable parameters, the A-type level 
scheme gives rise to self-phase modulation (SPM) due to 
the optical Kerr effect US]. As a consequence, the 
probe field experiences an intensity-dependent refractive 
index and thus acquires a time-dependent phase shift in 
accordance with its temporal intensity profile. This time- 
varying phase can be interpreted as a transient frequency 
shift. Furthermore, if the intensity of the probe pulse 
is sufficiently large, the SPM gives rise to characteristic 
modification of its Fourier spectrum. For example, an 
initially Gaussian Fourier spectrum typically break up 
into several peaks in the Fourier domain |43j . 
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A. Weak probe pulses 

In a first step we focus on weak probe pulses such that 
the effect of SPM is negligible. To this end, we eval¬ 
uate the Fourier spectra inside the medium and verify 
that typical manifestations of SPM are absent. In this 
parameter regime the remaining system parameters gov¬ 
erning the probe pulse propagation are given by A and 
Ac- In the following we discuss four different detuning 
conditions in order to systematically investigate the influ¬ 
ence of the Raman detuning A and Kerr detuning Ap/. 
For this, we evaluate all four combinations of A and Ap/ 


either 

0 or 4y: 





(*) 

A = 

0 

and 

Ap/ = 

= 4y, 

(ii) 

A = 

4y 

and 

Ap/ = 

= 4y, 

(Hi) 

A = 

0 

and 

Ap/ = 

= 0, 

(iv) 

A = 

4y 

and 

Ap/ = 

= 0. 


The resulting pulse propagation dynamics is shown in 
Fig.H It can be seen that Ap/ crucially affects the ab¬ 
sorption. The two cases with Ap/ = 0 lead to a trans¬ 
mission of less than 10% of the probe pulse intensity 
through the medium. In contrast, the cases Ap/ = dy 
give rise to approximately 75% probe intensity transmis¬ 
sion. This behavior can be traced back to the absorption 
on the |4) o |2) transition. Figure]^ shows the effect of 
A, which affects the shape and width of the transmitted 
pulse, whereas its effect on the transmitted intensity is 
small. 


B. Strong probe pulses 

Next, we consider probe pulses with higher intensities 
such that the effect of SPM is expected to become of 
relevance. Figure shows the probe pulse propagation 
dynamics as function of the propagation distance for case 
(ii), i.e., Ap/ = A = dy. The other parameters are cho¬ 
sen as in Figs. and except for the intensity higher 
by a factor of 10. Panel (a) shows the temporal probe 
pulse shape throughout the propagation. We find that 
the pulse gradually breaks up, until two fully separated 
pulses are obtained at r?C/7 = 100- The shapes of these 
sub-pulses remain Gaussian, and they are found to subse¬ 
quently propagate through the medium without further 
break-up. However, the individual subpulses broaden as 
they propagate through the medium, as it was also found 
in the corresponding weak probe pulse case o f Fig . 

To further analyze this pulse break-up, Fig. |4(b)] shows 
the power spectra of the probe pulse throughout its prop¬ 
agation, obtained by Fourier transformation of the tem¬ 
poral pulse shape. The power spectrum evolves from the 
Gaussian input spectrum into a double-peaked structure 
which is a characteristic signature of the SPM effect. 

A further numerical study reveals that the number of 
generated subpulses depends on the initial probe field 
intensity. Figure [^demonstrates that a three-peak pulse 


xlO’* 




(b) 


FIG. 4: (Color online) (a) Intensity |Gp/yp of a stronger 
probe pulse as function of the propagation distance. Param¬ 
eters are A = Ap/ = 4y and tip = \/0.15y; all other pa¬ 
rameters are chosen as in Fig.[^ (b) shows the corresponding 
probe pulse power spectra obtained by Fourier transformation 
as function of the propagation distance. 

can be generated from a single Gaussian probe pulse if 
the initial probe field intensity is suitably increased. As 
expected, the precise shape of the output pulse depends 
sensitively on the probe pulse intensity and the value of 
the detuning A. For example, the three-peak shape in 
Fig.j^turns into a two-peak structure if the three-photon 
detuning is increased to Api = A = 8y. Finally, we 
verified that the shaping and splitting of the probe pulse 
does not rely on the initial Gaussian pulse shape and can 
also be observed with non-Gaussian initial pulse shapes, 
such as secant-hyperbolic pulses. 

IV. THEORETICAL ANALYSIS 

Next, we provide a theoretical interpretation of the 
numerical results presented in Sec. |III[ To this end, we 
adapt the standard description of pulse propagation in 
optical fibers [33] to our system. 

We find that this model reproduces our numerical re¬ 
sults accurately for weak probe pulses where the SPM 
phenomenon is small. For parameters where the SPM ef¬ 
fect becomes dominant our model allows us to gain qual- 
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FIG. 5: (Color online) Propagation dynamics of a probe pulse 
with Up = Vo .57 and A — Api = 47 . All other parameters 
are chosen as in Fig. 


itative insights and provides an intuitive explanation for 


the pulse splitting phenomenon described in Sec. Ill 


Details on our simple model for the propagation of op¬ 
tical pulses in nonlinear media inspired by the standard 
fiber optics approach [43] are provided in Appendix 
There, we perform a perturbative analysis with respect to 
the probe field strength to derive the linear susceptibility 
x(a;) of the probe field, and find that the spatio-temporal 
evolution of the probe pulse can be described by 



FIG. 6: (Color online) (a) Real and imaginary parts of the 


linear susceptibility x in E<1- (A12l as a function of the 
Fourier component lu of the probe pulse envelope. We have 
Po = x(0) = 0, and j3n is the nth derivative of y at cj = 0, 
see Appendix The real and imaginary parts of P 2 [Pa] are 
shown in (b) (Jc)]. (d) Real and imaginary parts of R-p as a 
function of the Kerr detuning A / for A = 47 . In (a)-(d), 
blue solid (red dashed) lines show real (imaginary) parts and 
Dc = 7- 


pulse via 


_ 

dx~^ cdt 


^ VLp = ir] 


Vg = [ 1 /c-f 77 Re(/ 3 i)]’ 


(16) 




(14) 

The parameters Pn and Rp are determined by the micro¬ 
scopic model introduced in Sec. jllj In particular. 


Pn = 


d^X 

duj^ 


uj—0 


Ap—Ac 


(15) 


is the n-th derivative of the medium susceptibility x(w) 
evaluated at th e cen ter frequency w = 0 of the probe 
pulse [see Eq. (A13l]. Note that due to the approxi¬ 
mations used throughout the derivation, our model de¬ 
scribed by Eq. (14) is valid under the following condi¬ 
tions. (a) The probe pulse must be sufficiently long 
in time such that the atomic medium follows the probe 
pulse dynamics adiabatically. (/3) The probe field spec¬ 
trum must be sufficiently narrow such that the Taylor 
expansion of the linear susceptibility is justified [see Ap¬ 
pendix ^ for details]. ( 7 ) The probe field intensity must 
be weak enough such that all non-linear effects are well 
described by the third-order term proportional to Rp. 

We start by analyzing the relevant model parameters, 
see Fig. Panel (a) shows the linear susceptibility x{^) 
as a function of frequency uj scanned a cros s the probe 
pulse center frequency a; = 0 [see Eq. (A9)]. The pa¬ 
rameter Pq is the linear susceptibility x of the probe field 
evaluated at a; = 0, and vanishes for resonant Raman 
fields [see Eq. (HI)] if the ground-state decoherence rate 
r is neglected. Pi is related to the group velocity of the 


where Re denotes the real part. The parameters P 2 and 
P^ account for group velocity dispersion. The real and 
imaginary parts of P 2 [Ps] are shown in Fig. [^b) [(c)] as 
a function of the detuning A. Both real and imaginary 
parts of P 2 and P 3 determine the width and amplitude 
of the probe pulse at the output. For A = 0 only the 
imaginary part of P 2 is different from zero. This term 
accounts for absorption due to the finite width of the 
EIT window in frequency space. For larger values of A 
the real part of P 2 increases and describes the change of 
the group velocity for frequencies away from the central 
frequency of the pulse. In general, P 3 represents a higher- 
order correction to the group velocity dispersion. The 
nonlinearity due to the Kerr transition is proportional to 


Rp defined in Eq. (A 8 ), and the dependence of Rp on 


Ap/ is shown in Fig. [^d). The imaginary part of Rp 
is only large for Ap/ « 0 , and hence absorption on the 
Kerr transition is small for Apz/y ^ 1. Note that the 
quantitative impact of the term proportional to Rp in 


Eq. (14) depends on the probe pulse intensity. 


Equation (14) can be numerically solved using the 
split-step Fourier method [33|. We find excellent agree¬ 
ment with the numerical results obtained in Sec. IIIII for 
parameters where the model conditions (a)-( 7 ) are met. 
In particular, this is the case for all parameters where the 
SPM effect is small, and in this regime the explanation 
of all phenomena in terms of the parameters Pi and Rp 
is straightforward. We begin with a weak probe pulse 
and a large detuning Ap/ on the Kerr transition [cases 
(i) and (ii) in Sec. Ill . Since the imaginary part of Rp 
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FIG. 7: (Color online) Comparison between the full numerical 
solution to Eqs. (61 and (|^ [green solid line] and the simple 
model in Eq. ( 14 [[blue dot-dashed line]. The parameters in 
(a) are Up = 0.0257, Op = 50/7, A = 7 and Ap/ = 87.. 
Results in (b) correspond to flp = VO.I57, Cp = 30/7, s-nd 
A = Ap/ = 47. 


is small and hence the absorption on the Kerr transition 
is small for these parameters, the probe pulse can propa¬ 
gate through the medium without significant losses. The 
additional broadening of the probe pulse in case (ii) as 
compared to (i) can be explained by the non-zero group 
velocity dispersion for the parameters in (ii). Next we 
consider a resonant Kerr field Ap/ = 0 [cases (iii) and 
(iv) in Sec. Ill . The strong damping of the probe pulse 
found in Sec. Ill can be explained by the large imaginary 
part of Rp. The Kerr transition hence destroys the EIT 
phenomenon and leads to absorption of the the probe 
pulse. The non-vanishing group velocity dispersion in 
case (iv) leads to an additional broadening of the pulse 
as compared to (iii). 

In order to gain a qualitative insight into the pulse 
splitting phenomenon discussed in Secjlll B we identify 


a parameter regime where the pulse splitting occurs and 
the conditions (a,/ 3 , 7 ) of our model are met. This is 
shown in Fig. |7(a)[ where the pulse starts splitting at the 
end of the medium, and there is a reasonable agreement 
between the full numerical calculations and our theoret¬ 
ical model. From the model, the pulse splitting effect 
can be understood as follows. The SPM associated with 
Rp results in a splitting of the probe pulse in frequency 
space. In the presence of group velocity dispersion de¬ 
scribed by /?2 and /Sa, these pulses propagate at differ¬ 
ent speeds. This explains the emergence of two different 
pulses in the temporal domain. 

The pulse splitting as shown in Fig. is much more 
pronounced than in Fig. Our model predicts the pulse 
splitting effect for the parameters in Fig. but it is not 
sufficient for reproduc ing q uantitatively correct results. 
This is shown in Fig. |7(b)[ where we compare the the¬ 
oretical model to the full numerical simulation for pa¬ 
rameters as in Fig. |4(a)] The reason for the discrepancy 
between the simple model in Fq. (14) and the numerical 
approach is that the conditions ( 0 )^( 7 ) for the validity 
of the model are not met for the considered parameters. 
In particular, the probe pulse width in frequency space is 
not narrow enough to justify the adiabatic approximation 
rigorously. Second, the linear susceptibility y varies too 
strongly over the width of the pulse in frequency space 
for A = 47 such that the expansion of y up to third or¬ 
der is insufficient. Third, the probe pulse intensity is not 
small enough such that higher order terms beyond the 
third-order nonlinearity will contribute. 

Despite these shortcomings, we find that the qualita¬ 
tive explanation for the pulse splitting effect suggested 
by the simple model remains valid. To this end, we note 
that the power spectrum in Fig. Bb) found by numeri¬ 
cally fourier transforming the data in Fig. Qa) shows a 
clear signature of the SPM effect. The initially Gaus¬ 
sian pulse develops into two separate sub-pulses with in¬ 
creasing propagation distance. Let us denote the cen¬ 
tral frequencies of these pulses by uj+ and a;_, respec¬ 
tively. From the data in Fig. |^b) we find ui± « iO.OSy. 
Next we show that the splitting of the probe pulses 
in the temporal domain is associated with the different 
group velocities Vg{uj±) of these two sub-pulses. In or¬ 
der to verify this, we calculate Vg{u}±) from the disper¬ 
sion relation shown in Fig. according to Fq. ( [T^ 
with /3i evaluated at a; = uj±. As can be seen from 
Fig. Bb)i the two-pulse structure starts to develop at 
vC/l ~ 50 and is fully established at ryC/y ~ 100. We 
thus assume that the two sub-pulses evolve over a dis¬ 
tance TyA^/y « 130 with different group velocities un¬ 
til the end of the medium at jyC/y = 200 is reached. 
We estimate that the corresponding time difference is 
1307 / 77 |l/?;g(w_) — l/vg{uj+)\ « OO/y, which is in rea¬ 
sonable agreement with the time delay at ryC/y = 200 
according to Fig. ^a). In conclusion, we find that the 
pulse splitting arises from an interplay between SPM and 
group velocity dispersion. 

Finally, we note that the pulse splitting reduces the in- 
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tensity of the individual pulses, and suppresses nonlinear 
effects. Hence, the two pulses in Fig. 4(a)| do not split 
again. Higher intensities such as in Fig. lead to three 
pulses in the temporal domain. Interestingly, this result 
does not agree with the naive expectation that the pulse 
first breaks up in two pulses that subsequently split up 
in two pulses each, giving rise to a total of four pulses. 
This outcome is a manifestation of the non-linearity of 
the SPM phenomenon. 


V. CONCLUSION 

We have studied the propagation dynamics of a Gaus¬ 
sian pulse through an atomic medium in four-level N- 
type configuration under different input conditions of 
control and probe field intensities and detunings. As 
our first result, we have found that at two-photon res¬ 
onance condition with large Kerr detuning, distortion¬ 
less propagation of the probe pulse with low absorption 
is possible. In contrast, when the intensity of the probe 
pulse becomes higher, a splitting of the single Gaussian 
input pulse into a sequence of pulses may occur. To in¬ 
terpret the origin of this pulse splitting, we have derived 
a simple theoretical model for nonlinear pulse propaga¬ 
tion in our model system. The parameters of this model 
can be calculated from the system’s steady state den¬ 
sity matrix. We conclude that in particular an interplay 
between Kerr nonlinearity and group velocity dispersion 
leads to the pulse splitting. The nonlinearity leads to 
a decomposition of the probe pulse in frequency space. 
Subsequently, the such formed subpulses propagate with 
different velocity due to group velocity dispersion, such 
that a splitting in the time domain arises. The temporal 
shape of the probe pulse at the output does thus indi¬ 
cate the strength of the self-phase modulation induced by 
the Kerr-nonlinearity. Our numerical analysis has shown 
that the number of peaks can be controlled by adjusting 
the laser parameters. For example, the group velocity 
can be controlled by the detuning of the probe and con¬ 
trol fields with the excited state of the A sub-system. 
In addition, the Kerr-detuning and the laser intensities 
control the magnitude of the self-phase modulation. 
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Appendix A: Model parameters 


Here we outline the derivation of the model in Eq. (14) 
from Eqs. (§ and Q. We assume that the probe field is 
sufficiently weak (|r 2 p| <C 7 , |flc|) and pursue the pertur¬ 
bative approach outlined in [331 |3S] . The atomic density 


operator is expanded as 


00 

P = ^P^^\ (Al) 

where denotes the contribution to p in fcth order in 
the probe field Hamiltonian 


ifp = -(|2)(4|-|l)(3|)IIp + H.c.. 


(A2) 


In order to obtain the solutions p^^'> we re-write the mas¬ 
ter equation as 


^P ^qP ^ p] j 


(A3) 


where the linear super-operator Cq is independent of the 


probe field. Inserting the expansion (Al) into Eq. (A3) 


leads to the following set of coupled differential equations 




(A4) 


fc>0. (A5) 


Equation ( |A4[ ) describes the interaction of the atom with 
the control field to all orders and in the absence of the 
probe field. Higher-order contributions to p can be ob- 
tain ed if Eq. ( |A5[ ) is solved iteratively. Equations ( |A4[ ) 
and ( |A5[ ) must be solved under the constraints Tr(p(°iy = 
1 and Tr(p(*’^) =0 (k > 0 ). 

First we consider the nonlinear term where we employ 
the adiabatic approximation, i.e., we assume that the 
temporal length of the probe pulse is much longer than 
the lifetime of the excited states. In this case we can find 
the nonlinear response of the medium to the probe field 
by solving for the steady state of Eqs. (A4) and (A5) up 


to third order in Hp and obtain 


p(°) = |3)(3|. 


(A 6 ) 


Under the resonance condition A = Ap = Ac the only 
non-zero terms contributing in Eq. (j^ are P 24 p^i^. 
We find 


P^24 - Pl3 = . (A7) 

where 

2 

^" = "(Ap,+*7/2)|0,|2' 

Next we calculate the dominant correction arising from 
the spread of the probe pulse in frequency space. To 
this end we consider the Fourier expansion of the slowly 
varying probe pulse envelope. 


00 



— 00 











We solve Eqs. (A4) and (A5) up to first order and obtain 


The coefficients f3k are shown in Fig. and given by 


OO 

p[3it) = - j xM^pMe"*‘^*dw, (AlO) 

— OO 


where 

x(w) 


w + Ap — Ac 

{uj + Ap + ij/2){uj + Ap — Ac) — 

(All) 


o 

II 

o 

(A14) 


(A15) 

2(A + *7/2) 

(A16) 

, 6[|0,|2 + (A + *7/2)2] 

* “ |n,l» 

(A17) 


There are no other terms up to first order contributing 
to the coherences pi 3 or p 24 j and x is shown in Fig. 

If X is slowly varying over the typical width of the probe 
pulse, we can expand x ia a Taylor series around w = 0. 
For Ap = Ac = A and up to third order we obtain 


Inserting the expansion ( A12[) into Fq. (AIO) gives 


Pwit) 


g d’^Qp 

at" ■ 

n=0 


(AI8) 


where 


, ./£• 

x(w) « ^ —/ 3 „ , 
n\ 

n—O 


(A12) 


In summary, the right-hand side of Fq. Q can be written 
as 


iv{p24 - Piz) ~ ip 


(3) 

P2i 


(3) 

Pl3 


Pl 3 


(AI9) 


Pn 


d^X 


duj'^- 


uj —0 


Ar=Ac 


(A13) Subs tituti ng the results from Fqs. (|A7[ ) and (AI8) into 
Fq. (A19) allows us to obtain Fq. (Ildp. 


[1] M. Fleischhauer, A. Imamoglu, J.P. Marangos, Rev. 
Mod. Phys. 77 , 633 (2005). 

[2] S. E. Harris, Phys. Today 50(7), 36(1997). 

[3] L. V. Hau, S. E. Harris, Z. Dutton and C. H. Behroozi, 
Nature (London) 397, 594 (1999). 

[4] M. M. Kash, V. A. Sautenkov, A. S. Zibrov, L. Hollberg, 
G. R. Welch, M. D. Lukin, Y. Rostovtsev, E. S. Fry, and 
M. O. Scully, Phys. Rev. Lett. 82, 5229 (1999). 

[5] D. Budker, D. F. Kimball, S. M. Rochester, and V. V. 
Yashchuk, Phys. Rev. Lett. 83, 1767 (1999). 

[6] G. Heinze, C. Hubrich, and T. Halfmann, Phys. Rev. 
Lett. Ill, 033601 (2013). 

[7] K. P. Heeg, J. Haber, D. Schumacher, L. Bocklage, H.-C. 
Wille, K. S. Schulze, R. Loetzsch, 1. Uschmann, G. G. 
Panins, R. Riiffer, R. Rohlsberger, and J. Evers, Phys. 
Rev. Lett. 114 , 203601 (2015). 

[8] G. Liu, Z. Dutton, C. H. Behroozi, and L. V. Hau, Nature 
409 , 490 (2001). 

[9] D. F. Phillips, A. Fleischhauer, A. Mair, R. L. 
Walsworth, and M. D. Lukin, Phys. Rev. Lett. 86, 783 
( 2001 ). 

[10] T. N. Dey and G. S. Agarwal, Phys. Rev. A 67 , 033813 
(2003). 

[11] M. Fleischhauer and M. D. Lukin, Phys. Rev. Lett. 84, 
5094 (2000). 

[12] A. Andre and M. D. Lukin, Phys. Rev. Lett. 89, 143602 

( 2002 ). 

[13] M. Bajcsy, A. S. Zibrov, and M. D. Lukin, Nature 426, 
638 (2003). 

[14] F. E. Zimmer, A. Andre, M. D. Lukin, and M. Fleis¬ 
chhauer, Opt. Commun. 264 , 441 (2006). 

[15] F. E. Zimmer, J. Otterbach, R. G. Unanyan, B. W. Shore, 


and M. Fleischhauer, Phys. Rev. A 77 , 063823 (2008). 

[16] S. A. Moiseev and B. S. Ham, Phys. Rev. A 71 , 053802 
(2005). 

[17] S. A. Moiseev and B. S. Ham, Phys. Rev. A 73, 033812 
(2006). 

[18] Y.-W. Lin, W.-T. Liao, T. Peters, H.-G. Ghou, J.-S. 
Wang, H.-W. Cho, P.-C. Kuan, and 1. A. Yu, Phys. Rev. 
Lett. 102 , 213601 (2009). 

[19] G. Nikoghosyan and M. Fleischhauer, Phys. Rev. A 80 , 
013818 (2009). 

[20] M. Fleischhauer, J. Otterbach, and R. G. Unanyan, Phys. 
Rev. Lett. 101 , 163601 (2008). 

[21] M. Kiffner and M. J. Hartmann, Phys. Rev. A 82 , 033813 

( 2010 ). 

]22] S. E.Harris, J. E. Field, A. Imamoglu, Phys. Rev. Lett. 
64, 1107 (1990). 

[23] H. Wang, D. Goorskey and M. Xiao, Phys. Rev. Lett. 87, 
073601 (2001). 

[24] J. P. Marangos, J. Mod. Opt. 45, 471 (1998). 

[25] R. Fleischhaker and J. Evers, Phys. Rev. A 77 , 043805 
(2008). 

[26] H. Schmidt and A. Imamoglu, Opt. Lett. 21, 1936 (1996). 

[27] H. Kang and Y. Zhu, Phys. Rev. Lett. 91, 093601 (2003). 

[28] G. F. Sinclair and N. Korolkova, Phys. Rev. A 76, 033803 
(2007). 

[29] J. H. Li, X. Y. Lu, J. M. Luo, and Q. J. Huang, Phys. 
Rev. A 74 , 035801 (2006). 

[30] J. Sheng, X. Yang, U. Khadka and M. Xiao, Opt. Exp. 
19 , 17059 (2011). 

[31] M. Bajcsy, S. Hofferberth, V. Balic, T. Peyronel, M. 
Hafezi, A. S. Zibrov, V. Vuletic, and M. D. Lukin, Phys. 
Rev. Lett. 102 , 203902 (2009). 




















9 


[32] M. Kiffner and M. J. Hartmann, Phys. Rev. A 81 , 
021806(R) (2010). 

[33] M. Hafezi, D. E. Chang, V. Gritsev, E. A. Dernier and 
M. D. Tonkin, EPL 94, 54006 (2011). 

[34] D. E. Chang, V. Gritsev, G. Morigi, V. Vnletic, M. D. 
Lukin, and E. A. Dernier, Nat. Phys. 4, 884 (2008). 

[35] A. Andre, M. Bajcsy, A. S. Zibrov, and M. D. Lukin, 
Phys. Rev. Lett. 94, 063902 (2005). 

[36] L. E. Zohravi, M. Abedi and M. Mahmoudi, Commun. 
Theor. Phys. 61 , 506 (2014). 

[37] T. N. Dey and G. S. Agarwal, Phys. Rev. A 76, 015802 
(2007). 

[38] S. E. Harris and Y. Yamamoto, Phys. Rev. Lett. 81 , 3611 
(1998). 

[39] D. Maxwell, D. J. Szwer, D. Paredes-Barato, H. Busche, 
J. D. Pritchard, A. Gauguet, K. J. Weatherill, M. P. A. 


Jones, and G. S. Adams, Phys. Rev. Lett. 110, 103001 
(2013). 

[40] S. Sevincli, N. Henkel, G. Ates, and T. Pohl, Phys. Rev. 
Lett. 107 , 153001 (2011). 

[41] C. S. Hofmann, G. Gunter, H. Schempp, M. Robert-de- 
Saint-Vincent, M. Garttner, J. Evers, S. Whitlock, and 
M. Weidemiiller Phys. Rev. Lett. 110 , 203601 (2013). 

[42] A. V. Gorshkov, J. Otterbach, M. Fleischhauer, T. Pohl, 
and M. D. Lukin, Phys. Rev. Lett. 107 , 133602 (2011). 

[43] G. P. Agrawal, Nonlinear Fiber Optics (Academic, 
Boston, 1989). 

[44] M. Kiffner and K.-P. Marzlin, Phys. Rev. A 71 , 033811 
(2005). 

[45] M. Kiffner and U. Dorner and D. Jaksch, Phys. Rev. A 
85 023812 (2012). 



